Efficient extraction of quantum Hamiltonians from optimal laboratory data 
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Optimal Identification (01) is a recently developed procedure for extracting optimal information 
about quantum Hamiltonians from experimental data using shaped control fields to drive the system 
in such a manner that dynamical measurements provide maximal information about its Hamiltonian. 
However, while optimal, OI is computationally expensive as initially presented. Here, we describe 
the unification of 01 with highly efficient global, nonlinear map- facilitated data inversion procedures. 
This combination is expected to make OI techniques more suitable for laboratory implementation. 
A simulation of map-facilitated OI is performed demonstrating that the input-output maps can 
greatly accelerate the inversion process. 

PACS numbers: 42.55.-f 



I. INTRODUCTION 

A general goal in atomic and molecular physics is to 
quantitatively predict quantum dynamics from knowl- 
edge of the system Hamiltonian. However, sufficiently 
accurate information about these Hamiltonians is still 
lacking for many applications. Although the capabilities 
of ab initio methods are improving, they remain unable 
to provide the quantitative accuracy needed to predict 
many quantum dynamical phenomena and data inver- 
sion remains the most reliable source of precision infor- 
mation about quantum Hamiltonians. But, traditional 
data inversion techniques are hindered by the fact that 
(a) spectroscopic and collision data only provide informa- 
tion about limited portions of the the desired interactions 
and (b) the relationships between quantum Hamiltonians 
and their corresponding observables are generally nonlin- 
ear pj. 

The recently proposed optimal identification (OI) 0,0 
procedure provides a new approach to quantum system 
Hamiltonian identification through laboratory data in- 
version. The operating principle behind OI is to improve 
the information content of the data by driving the quan- 
tum system using a tailored control field (e.g., a shaped 
laser pulse). If suitably chosen, the control field forces the 
data to become highly sensitive to otherwise inaccessible 
portions of the Hamiltonian and therefore enables a high 
fidelity inversion. A key component of the OI concept ex- 
ploits the fact that the inversion will generally produce a 
family, or distribution, of Hamiltonians that are consis- 
tent with the laboratory data 0, Q . The breadth of this 
inversion family provides a figure of merit for the quality 
of the inversion. The limiting experimental factors that 
prevent the inverted family of Hamiltonians from collaps- 
ing down to a single (i.e., completely certain) member 



arise from two sources. First, the finite precision of the 
data reduces the resolving power of the measurements 
and makes it possible for multiple Hamiltonians to be 
consistent with the data to within its experimental error. 
Second, Hamiltonians that differ in ways for which the 
data is insensitive also reduces the inversion quality. OI 
operates by attempting to drive the quantum system into 
a dynamical state where the associated experimental er- 
rors are least compromising, yet where the measurements 
provide distinguishing power between Hamiltonians to 
narrow down the family consistent with the data. 



As originally formulated, the OI algorithm requires 
performing a large number of global, nonlinear data in- 
versions, each of which can be computationally expen- 
sive Q- Resolving the members of the inversion family 
consistent with the data might require extracting hun- 
dreds of distinct Hamiltonians involving numerous so- 
lutions of the Schrodinger equation. In this paper, we 
demonstrate that it is possible to greatly reduce the com- 
putational challenge of performing OI by incorporating 
map-facilitated inversion techniques^, 0, E| that have 
been specifically developed for efficiently finding these 
solution families. A map is a predetermined quantita- 
tive input^output relationship which can alleviate the 
expense of repeatedly solving the Schrodinger equation 
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This paper provides a detailed description of the al- 
gorithm demonstrated in Ref. for extracting both in- 
ternal Hamiltonian and transition dipole moment matrix 
elements from simulated laser pulse shaping and popula- 
tion data. Section UTI reviews the OI concept introduced 
in Ref. Q and then extends this procedure to incorporate 
map- facilitated inversion. Section IlIII provides a detailed 
description and in-depth analysis of the simulations that 
were presented in Ref. 0- 
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II. ALGORITHM 

The principal behind the 01 algorithm is to optimize 
an external control field that drives the system in a 
manner similar to the learning-loop techniques utilized 
in many current coherent quantum control experiments. 
The distinction between 01 and the latter control ex- 
periments lies in the optimization target: 01 is guided 
to optimize the extracted Hamiltonian information. The 
01 learning algorithm programs a control pulse-shaper 
to drive the quantum system followed by the collection 
of associated dynamical observations. In practice, the 
control optimization is performed over a discrete set of 
variables (control "knobs"), c = {ci, . . . ,Cat c }, 



E(t) -> E(t;ci,c 2 , ■ ■ ■ ,c N J 



(1) 



where the space of accessible fields is defined by vary- 

n (vain) ~ ~ (max) 

mg each Ci over a range , c\ < c, < c\ 



For each trial field, Ek(t;ck), k — 1,2,..., a set 
of measurements, <& k lab \ are performed on the sys- 
tem. Each trial field yields M individual measure- 
ments (e.g., the populations of M different quantum 
states), 4>^ ab ) = {$ k l< f\ . . . , ' }' w ith associated er- 
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Inversion is performed by adopting a discrete set of 
variables, h = {hi, ■ . ■ , fiN h } used to distinguish one 
trial Hamiltonian from another. There are many pos- 
sible ways to define these Hamiltonian variables, and the 
best representation must be selected to suit the quan- 
tum system being inverted, but in general, a sufficiently 
flexible and accurate description of Hamiltonian space re- 
quires a large number of variables, N n 3> 1. Inversion is 
accomplished by minimizing, 
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where $ fc . m [h] is the m th (m = 1,...,M) observable's AH*[E k (t)], 
computed value for the trial Hamiltonian, h, under the 
influence of the external field, Ek(t). Optionally, a regu- 
larization operator, K, acting on the Hamiltonian, h, can 
be used to incorporate a priori behavior, such as smooth- 
ness, proper asymptotic behavior, symmetry, etc., into 
the inverted Hamiltonian 0, Q, 13 ■ While the data error 
distributions are assumed to have hard bounds, in 
Eq. J2J), other distributions could be used as well. 

The output of the inversion optimization is a set of N s 
Hamiltonians, {h*, . . . , }, that each ideally reproduce 

the measured observable, (j>^ afc ^ ) to within its experimen- 
tal error. The upper and lower bounds of each inverted 
variable defines the family, 



<h* = min{/i*J 



= max{/i*J 



(3) 
(4) 
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AH*[E k (t)} = _]Tj inv [h:;$f b) ] 
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where h* is the s th member of the inversion family found 
from Ek(t) and J7 mv is given by Eq. (J2J). The first term 
in Eq. © measures the ability of the inversion family to 
reproduce the data and the second measures the inversion 
uncertainty with a > being a coefficient that balances 
them. 

This measure of inversion uncertainty is used to guide 
the control optimization where the objective is to opti- 
mize AH* over the space of accessible fields by minimiz- 
ing the control cost function, 



where h* i is the I 



Hamiltonian variable from the s 



ih 



member of H*. The uncertainty in each Hamiltonian 
variable, Ah*, is quantified by the width of its corre- 
sponding solution space, 



Ah* =>h1 



(5) 



and the width of the family for each Hamiltonian variable 
is used to compute the uncertainty in the full inversion, 



J c [E(t;c)} = AH*[E(t;c)} 



N r 



Ci - c 
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(7) 



where the first term reduces inversion error and the sec- 
ond removes extraneous field components relative to their 
minimum value, balanced by (3 > 0. 

The result of the control field optimization is a set of 
laboratory data and its corresponding inversion results, 
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{El(t), $£, e* k , h* k } that provides the best possible knowl- 
edge of the unknown Hamiltonian provided by the set of 
accessible optimal control fields. This inversion result 
provides the Optimal Identification of the quantum sys- 
tem with uncertainty, 

Ahi*= >h* -<h* (8) 

where 'Ti,* and are computed using Eqs. J3) and Q 
for the optimal data. 

The most expensive operation in 01 is the many solu- 
tions of Schrodinger's equation required in the inversion 
component of the algorithm. The map-facilitated inver- 
sion aims to alleviate this costly operation. Computing 
quantum observables, 3?fe[h], from a given Hamiltonian, 
h, and control field, Ek(t; c^), defines a forward map, 

/ : (h;£ fc )-**M (9) 

that is parameterized by the control field, Ek(t). In the 
initial presentation of 01, this map is explicitly evaluated 
each time a trial Hamiltonian is tested to determine if it 
is consistent with the laboratory data. However, it has 
recently been found that it is possible to pre-compute 
this map to high accuracy by sampling f(h;Ek) for a 
representative collection of Hamiltonians. 

Although it is generally impossible to resolve / on a full 
grid in h-space due to exponential sampling complexity 
in Nh 3> 1 dimensions, it has been found [l| that an 
accurate nonlinear map can often be constructed using 
the functional form, 

N h N h 

/(h; E k ) = f (E k ) + fi(M; E k ) + fa(hi, h j; E k ) 

i—1 i<Cj 

+ ... + f 1 ... Nh (h 1 ,...,h Nh ;E k ) (10) 

where fo(E k ) is a constant term, the functions, 
fiihi) E k ), dependent upon a single variable, hi, the func- 
tions fij(hi, hj; Ek) depend upon two variables, hi and 
hj, etc. Expansions of the form in Eq. IjlOjl belong to a 
family of multivariate representations used to capture the 
input— >output relationships of many high-dimensional 
physical systems 0, U [3 (3 13 111 13 111 G3 It has 
been shown that Eq. IjlOjl converges to low order, L <C Nh 
for many Hamiltonian— s-observable maps. A low order, 
converged map expansion can be truncated after its last 
significant order without sacrificing accuracy or nonlin- 
earity which dramatically reduces the computational la- 
bor of constructing the map. The specifics for how the 
individual expansion functions should be evaluated can 
be found in previous papersp],H^|. 

Map-facilitated 01 operates by replacing the step of 
explicitly computing $fc[h] in Eq. with the value 
f(h;Ek) obtained by evaluating the functional map in 
Eq. IjlOjl which serves as a high-speed replacement for 
solving the Schrodinger equation. The map (or poten- 
tially multiple maps if this is necessary to ensure suffi- 
cient accuracy |4j) must be constructed as a preliminary 



step prior to initiating the inversion optimization. Since 
map construction requires knowledge of the control field, 
it may not be possible to pre-compute all of the maps 
prior to executing the full 01 loop. Instead, a new map 
will generally need to be constructed for each trial laser 
pulse. 

III. ILLUSTRATION 

The map-facilitated 01 algorithm was simulated for an 
8-level HamiltonianpjJ chosen to resemble vibrational 
transitions in a molecular system where the objective 
was to extract optimal information about the molecu- 
lar Hamiltonian, H, and dipole moment, fi, for a system 
having the total Hamiltonian, 

H = H-fi-E k (t). (11) 

The control fields have the form, 

25« (t) = exp ( ~ ( ^J /2)2 ) £ A ? } cos (^ + °i ] ) 

(12) 

where the uji are the resonance frequencies^^ °f A\ f , 
their corresponding amplitudes and Of , their associated 
phases. Control field noise was modeled as parametric 
uncertainties in the Ai and 

A^ = (l + 7 ^)A h ^ = (1+7^ (13) 

where different random values between ±e^ id ) were cho- 

(i) (i) 
sen for 7^ and for each pulse. 

The 01 simulation involved learning the matrix ele- 
ments of the Hamiltonian, H pq — (p\H\q), and dipole mo- 
ment, \i vq = (p\fi\q) for a chosen basis, \p), p = 1, . . . , 8 
|l8j. Simulated laboratory data was generated by prop- 
agating the initial wavefunction under the influence of 
the applied control field from its initial state and com- 
puting the populations (in the chosen basis, \p)) at var- 
ious times, t q , during the evolution. Each population 
"measurement" was averaged over D = 100 replicate ob- 
servations for a collection of noise-contaminated fields, 
{E k \t)}, j — 1,...,D, centered around the nominal 
field, Ek(t) to simulate experimental uncertainty in the 
laser pulse-shaping process. Measurement error, e ( - lab ) 1 
was introduced into the population observations accord- 
ing to, 

<S=(^ + P^Pi(t q -AE^ (*)])). =i d (14) 

where the pij were chosen randomly between ±e^ obs \ the 
relative error in each population observation. A different 
random value, pij, was selected for every simulated mea- 
surement. 

Equation Q was minimized over c = {Ai, 6{\ using a 
steady state GA with a population size of 30, a mutation 
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FIG. 1: (color online) Comparison of optimal versus conventional Hamiltonian identification. In plots (A-C) the identification 
error (darker shading implies larger inversion error) respectively reflects OFs performed using 1, 2, and 4 samples of the 
populations during the wavepacket evolution. Increasing the number of data points available to the OI process improves the 
quality of the extracted Hamiltonian. Plot (D) represents a conventional identification where the populations were sampled 25 
times during the time evolution. Despite using significantly less data, OI obtained higher quality Hamiltonian information. 



rate of 5%, and a cross-over rate of 75%. The pulse 
parameters where chosen to be T — 1.0 ps and s = 200 fs, 
the amplitudes, Ai, were allowed to vary over [0,1] V/A, 
and the phases were allowed to vary over [0,27r] rad. The 
laboratory measurement error was assumed to be £^ a ^ = 
2%, the field error, e^ ld "> = 1% and the observation times, 
t q , were uniformly spaced over the evolution period (i.e., 
the time between observations was At = T/Q with Q — 
1, 2 or 4). At each time, the full set of 8 populations was 
measured. The parameters, a in Eq. and (3 in Eq. 
0, were ramped from lxlO" 4 to lxlO~ 2 over the GA 
evolution, although the optimization was insensitive to 
the exact choices of a and (3. Typically ~50 generations, 
or approximately 800 trial fields, were needed for GA 
convergence. 

Global inversion to identify the Hamiltonian family 
corresponding to the data, $^ ab \ was performed by min- 
imizing Eq. J2)l using a map-facilitated inversion algo- 
rithm. For each inversion, the family of consistent Hamil- 
tonians was identified using a steady-state GA with a 
cross-over rate of 70% and a mutation rate of 5%. The 
trial family size was N s = 500 and the GA population size 



was N p — 100. The Hamiltonian-space map variables, 
hi, were the matrix elements of the molecular Hamilto- 
nian, H mn , and the the dipole, /i mn . For the eight-level 
system, there were 36 Hamiltonian elements (symmetric, 
upper-triangle including the diagonal) and 28 transition 
dipole moments (symmetric, upper-triangle without the 
diagonal) producing an N% = 64-dimensional map. All 
maps were constructed to first order, L = 1, and 5* = 6 
sample points were used to resolve each map function for 
interpolation. The Hamiltonian-space domain extended 
±30% around its nominal value (each matrix element was 
assumed known to ±15% prior to the present identifica- 
tion). Typical map construction required an average of 
84 seconds to perform on an SGI MIPS single processor 
machine, and map-facilitated inversion required an av- 
erage of 51 seconds to converge. A single evaluation of 
the Hamiltonian— ^population map typically required ~ 1 
ms, while a similar solution of the Schroinger equation for 
this system took ~ 2 s. This difference is the origin of 
the savings associated with map-facilitated OI. 

The performance of the map-facilitated OI algorithm 
was assessed with the following four tests: 
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(A) An 01 was performed using populations measured 
at Q = 1 time, ti = T, producing 8 observations 
for the 64 unknowns. 

(B) An 01 was performed using populations measured 
at Q = 2 times, t\ = T/2 and t 2 — T, producing 
16 observation for the 64 unknowns. 

(C) An 01 was performed using populations measured 
at Q — 4 times, t q — qT/4, q — 1, . . . , 4, producing 
32 observations for the 64 unknowns. 

(D) A conventional inversion (with a randomly selected 
field) was performed using populations sampled at 
Q = 25 times, t q = qT/25, q = 1, ... ,25, producing 
200 observations for the 64 unknowns. 

The power spectra of the optimal control fields for the 
Q = 1, 2, and 4 01 inversions are shown in Figure Q] 
along with a graphical depiction of the 01 error, A/ij*. 
For a single time sample, Q = 1, the overall inversion 
error, computed as the average, 2(Ahi* + > hj* 
)), was found to be 3.67%. Figure E^A) shows that the 
majority of the error was contained in and around the 
diagonal elements of H (note that the Hamiltonian is not 
diagonal in the chosen basis, \p)). The average relative 
error in the dipole, 0.9413%, was significantly smaller, 
and the majority of dipole uncertainty appeared in the 
v— >v + 2 elements. The inversion error for the two-point, 
Q = 2, 01 demonstration was reduced to 2.910% for 
the molecular Hamiltonian. Again, the majority of the 
inversion uncertainty resides in the diagonal elements of 
H. The transition dipole moment elements also improved 
with an overall average uncertainty of 0.6270%. 

The inversion error, A/ij*, for a simulated 01 utiliz- 
ing data at Q = 4 times, is essentially eliminated. The 
average uncertainty in both the molecular Hamiltonian 
elements and in the transition dipoles is an order of mag- 
nitude smaller than the simulated error in the data, el' °^ 
= 2 %. The most dramatic demonstration of the OI's ca- 
pabilities is seen by comparing Figure and D). The 



plot in (D) represents a conventional inversion, performed 
using a randomly selected field and Q — 25 time samples, 
compared to Q = 4 for 01. The conventional inversion 
therefore had access to M — 200 data points while the 
map- facilitate 01 demonstration only had M = 32. De- 
spite what would appear to be a significant advantage in 
the amount of available data, the conventional inversion 
displays greater than two orders of magnitude more er- 
ror. The conventional dipole moment inversion produced 
an average uncertainty of 19.5% and the precision in the 
molecular Hamiltonian was 28.1%. The map- facilitated 
01 found a control field and associated data that essen- 
tially prevented the laboratory noise from propagating 
into the identified Hamiltonian information. 

IV. CONCLUSION 

We have presented simulated experimental data that 
demonstrates the utility of map-facilitate Optimal Iden- 
tification for the real-time laboratory identification of 
quantum Hamiltonians from dynamical physical observ- 
able data. The central concept behind 01 is that by suit- 
ably driving the quantum system in an optimal manner, 
it is possible to extract high-precision information about 
all relevant aspects of its Hamiltonian despite finite lab- 
oratory error in both the external control fields and mea- 
sured data. In this demonstration of map-facilitated 01, 
it was possible to improve the computational efficiency 
of data inversion by more than an order of magnitude 
compared to the initial presentation of the 01 concept 
0- This increased efficiency is ultimately expected to 
aid, if not be essential, for the practical implementation 
of 01. 
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